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Abstract 

We present unitarity as a method for determining the infinities present in graviton 
scattering amplitudes. B The infinities are a combination of IR and UV. By understanding 
the soft singularities we may extract the UV infinities and relate these to counter-terms 
in the effective action. As an demonstration of this method we rederive the UV infinities 
present at one-loop when gravity is coupled to matter. 
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1. Introduction 



Calculations in perturbative gravity [1] are notoriously difficult to perform. In par- 
ticular determining the renormalisability of gravity whether coupled to matter or not is 
a difficult issue. General arguments regarding the symmetries in the action may allow or 
prohibit counter terms. Such arguments show that the UV infinities vanish on-shell up 
to one-loop for pure gravity amplitudes and for up to two-loops for the particular matter 
coupling in supergravities [2]. However in gravity coupled to general matter, there are 
possible counter terms at one-loop level and a calculation must be performed to determine 
the coefficient of the potential counter-term. Such calculations have been done for gravity 
coupled to matter and the explicit coefficients obtained [3,4,2]. These calculations use an 
algorithm due to 't Hooft and Veltman [3] which examine counter terms in the effective 
action using the background field method [5]. (The coefficients of the counter terms in 
pure gravity at two loops have been calculated [6,7] with non-zero result.) 

Recently developments have been made in the calculation of on-shcU amplitudes using 
string inspired techniques and the Cutkosky rules. These have enabled new results to be 
obtained both in QCD [8] and perturbative quantum gravity [9,10]. In this paper we show 
how these techniques may be used to detect the effect of the counter terms in the one- 
loop amplitudes - and in fact extract their coefficient. This calculational method is quite 
distinct from usual techniques. 

The Cutkosky rules [11] relate two amplitudes 'sewn' together to the discontinuous 
parts of an amplitude at higher-order. Knowing the lower order expressions we may impose 
constraints upon the higher order parts of the amplitude. (For example, the Cutkosky rules 
relate the imaginary parts of one-loop amplitudes to tree amplitudes sewn together.) The 
constraints these rules impose upon amplitude can completely reconstruct some amplitudes 
and constrain the expressions for others [12,13]. In gravity this technique has provided a 
check for the 1-loop results from the string-inspired rules although several of the four point 
amplitudes could have been reconstructed entirely from the cuts. In general it is the "more 
supersymmetric" amplitudes which are cut-constructable. For general amplitudes we may 
construct, via the Cutkosky rules, expressions containing the correct cuts in all channels. 
Although such expressions contain finite ambiguities, we may use these expressions to 
evaluate the infinities present in an amplitude. In fact, we will demonstrate the one- 
loop UV-infinities present in matter-coupled gravity using this method. To enable us to 
disentangle the UV and IR infinities we also determine the general structure of the one-loop 
soft infinities. 

Finally, although amplitudes naively calculated using the Cutkosky rules contain finite 
ambiguities these may be resolved by a more sophisticated use of the rules [14]. We 
illustrate this with a specific calculation at one-loop involving graviton scattering in a 
theory with massive scalar matter. 
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2. Infrared divergences 



As we will see later, one can, by using unitarity, extract the infinities present in 
an amplitude. However, this method will not distinguish between UV and IR infinities. 
Nonetheless, by knowing the expected form of the IR singularities we may identify the 
remaining UV infinities. ^ 

In this section we determine the IR singularities present in amplitudes involving mass- 
less particles. The IR singularities in a QCD amplitude have been calculated [15] and our 
determination of the IR singularities for gravity follow in a similar manner. One may 
also deduce the form of the infra-red divergences using general arguments to imply uni- 
versality and then extract the form from a specific amplitude [16]. As is well known [17] 
IR singularities in a on-shell loop amplitude occur in soft limits of the loop momentum 
integrals. That is, when the momentum of internal propagators goes to zero. However, 
not all such occurances produce IR singularities. Any masslcss propagator through which 
loop momentum flows will be singular for a speciflc value of the loop momentum but this 
does not usually yield a singularity. To see this note that, in four dimensions, the singular 
part of the momentum integral can be expressed 



which is finite at the lower range of integration, \p\ — > the singularity having been sup- 
pressed by the factor of \p\^. (We always choose gauges where the propagator is Feynman- 
like ~ 1/p^. Gauges closely related to the string based rules [18] and the "World-line" 
approach [19] have this feature.) In general three adjacent propagators must vanish simul- 
taneously to obtain a soft divergence. To see this note that three adjacent propagators in 
a loop will be of the form (we shift the loop momenta such that the middle momenta is 
the integration momenta) 



We may obtain an IR singularity provided = and = 0. One can imagine more 
complicated situations in multi-loop diagrams where adjacent propagators vanish at choices 
of the multiple integrations. However the equivalent factors to in the above equation 
are of higher powers and suppress the infinities more strongly. After inspection of the 
possibilities we find the above situation is the generic case for a soft divergence in the 
loop-momentum integral. 

We will only, for general kinematics, find propagators where = and iiTf = if 
Ki and K2 are adjacent external on-shell momentum (for massless particles.) The middle 

^ We thank Lance Dixon [16] for help in reahsing this possabihty 




(2.1) 



111 1 1 1 



(2.2) 



{p + KiY p^ {p - K2Y (p2 + 2j9 • + Kl) p2 (p2 _2p.K2 + Kl) 
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propagator is then joining together the two external legs. One may see this general 
arrangement in fig. 1. Thus the soft limit may be found by taking the amplitude with 
one less loop and adding a soft particle between two of the external lines. Consider first 
an individual diagram, rather than the whole amplitude. We initially restrict ourselves 
to a diagram containing only gravitons. Throughout, amplitudes will be calculated using 
dimensional reduction [20]. 

Consider adding a soft graviton to an graviton scattering diagram by attaching a 
graviton between external legs with momentum ki and /c2 as in fig. 1. In doing so we must 
add a propagator for the graviton, two three point vertices and two extra propagators 
with momenta ki -\- p and /c2 — P to the normal Feynman diagram expression. The soft 
singularity, as in eqn. (2.2), will occur in that region of the J d^p integral where p is close 
to zero. 



kj+p 




k2 kj- p 

Figure 1. Representation of the soft limit calculation. 

Initially consider the soft-divergences in a one-loop diagram. Thus the soft graviton 
leg will connect two external legs of a tree diagram. Take the off-shell tree diagram as 
a function of the two external momentum ki and /c2 only and of the indices on these 
two legs only (keeping all others fixed), ^(|^^^) (-y^)(^i5^2) The on-shell tree diagram is 
the contraction of this with the external polarisation tensors (together with ki and /c2 
becoming onshell). 



6i"^e2-^^A;-^\(^^)(fci,fc2) (2.3) 

the amplitude being the sum of all such diagrams. The 1— loop diagram can be written in 
terms of the off-shell tree diagram with the additional three propagators and two vertices, 

^4-2e^ T("^)'('^^)(A;i,/C2,p) 



A 



tree 



,Aki+p, k2-p) 



(27r)4-2^ {p -t- /ti) V(P - ^ ("^)'(^'^) 
The tensor T is 

T("«'('^-)(fcl, k2,p) = ei"^^^e2^^"^ V(,,/3,),(57),(a/?)(^2, -p, "^2 + p) 

X P{6-y),{ap) V(^f,^^^)^^ap),{t,u)ikl:P:-kl - p) 



(2.4) 



(2.5) 
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where V(juz/),(CTp),(7r;) is the 3-graviton vertex and P(af3),{ap) is the propagator [1,21,7]. In 
the p — > hmit, we find T becomes 

T(«/3)'(7p)(/,^^/,2,p) = -Y^i"^e2^''(A;i +A;2)^ + 0(p) (2.6) 
To examine the leading soft singularity we must therefore look at the integral, 

^ (2.7) 



(2.8) 



(27r)4-2^ {p + kiYp^{p-k2Y 
which can be evaluated by usual Feynman parameter methods with the result, 

— zrr 1 



(47r)2-6(_s^2)i+6 e 
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where rr = r^(l — e)r(l + e)/r(2 — e) and Si2 = (/ci + /c2)^- ^ The leading singularity as 
p — > is then simply a factor multiplying the tree diagram. Summing over all diagrams 
then gives a factor multiplying the tree amplitude. To be precise, we find that the IR 
divergence, due to a soft graviton exchange between external legs 1 and 2, in a one-loop 
graviton scattering amplitude is 

(5^^ X (-.i^)-' X (2.9) 

There are a few subtleties with the above result. Firstly, we have to take the limit of 
an off-shell amplitude A^^^^{ki + p, k2 — p) as p ^ 0. This is trivial for a tree amplitude 
but for the n-loop case is more subtle. In general this will yield the on-shell result for an 
amplitude although this is not true merely for a single n-loop diagram [24] . Secondly there 
are sub-leading terms in eqn.(2.6). These vanish on summing over diagrams because we 
are then dealing with a physical on-shell amplitude. 

We have calculated the IR divergence due to soft graviton exchange between a specific 
pair of external legs. The total IR divergence in a one- loop amplitude is simply the above 
result summed over all pairs of external legs. 

2 ^ 

^rr ^ X ^^-^^^ X V(-Si,)i-^ (2.10) 



(47r)2-^ 4e2 ^ 

In gauge theories, amplitudes contain IR divergences from the self-energy corrections to 
the external legs. [15]. However, by power counting, such divergences do not occur for 
external graviton legs. 



^ We shall always give amplitudes (and integrals) in the unphysical region where all momentum in- 
variants such as Sij are negative. One may continue back to the physical regime by ln(— s) 
ln(|s|) — iTT@{s) etc. We also use the convention that all particles are outgoing. 
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Consider a specific example: namely that of a 1-loop four graviton amplitude. These 
are given in ref. [10]. Also let us look at the specific helicity configuration ^^~^°°p(1~, 2~, 3+, 4+). 
Summing over all pairs of legs the expected IR divergence is 



2{i-sy-^ + i-ty-^ + i-uy-^) x ^4^^^'"^(l-,2-,3+,4+) (2.11) 



(A factor of 2 arises because Si2 — S34 = s.) This is 



irr^'^ ( {s ln(— s) + t ln(— t) + u ln(— w)) 



^tree^^-^2-,3+,4+) (2.12) 



For this amplitude we expect no UV singularities from the general result that pure gravity 
is one-loop UV finite. If we examine the complete result [10], 



^i-ioop(i-^ 2-,3+, 4+) = ^^^^A--(l-, 2- 3+, 4+) 



X 



2/ln(-w) In(-t) ln(-s)\ ^ . \ 

- ^ + ^ — + ^ — + finite terms 

e\st su tu J J 



(2.13) 



we find the expected IR infinity structure. 

In amplitudes with scalars and gravitons, there are also soft divergences which can 
be calculated similarly. For the case of a scalars and gravitons there are six configurations 
which must be considered, as shown in fig. 2. In this figure scalars lines are dashed lines 
whereas graviton lines are solid. The analysis in all cases closely follows the pure gravity 
calculation. We find that (a), (b) and (c) all give the same result as the pure gravity case, 
(2.10); (d) and (e) give no soft contribution. The case (a), (b) and (c) are the diagrams 
where a soft graviton is exchanged between the two external legs. The result in eqn. (2.10) 
is then universal. That is, true whether external legs are gravitons or scalars. We can 
also consider the IR divergences with other types of external particle and find the same 
universal nature of the soft divergences. This universality may be expected upon general 
arguments [16]. 
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3. Divergences from the Cutkosky Rules 

As is well known, [11], unitarity in the form of the Cutkosky rules, relates the imagi- 
nary part of a one-loop amplitude to the product of two tree amplitudes integrated over all 
intermediate states. In ref.[13] situations where the Cutkosky rules determine the ampli- 
tude entirely (and not merely the imaginary part ) were investigated with the result that 
in supersymmetric gauge theories the amplitudes may be determined by the Cutkosky 
rules alone. For gravity only a few amplitudes may be determined from unitarity alone [8] . 
However, we may determine the infinity structure entirely from the Cutkosky rules. 

In general, we consider the cut in the channel (/cmi + ^mi+i + • • • + km2-'i. + km^)'^ for 
the 1-loop amplitude A{1, 2, . . . , n), depicted in fig. 3 and given by 

J rfLIPS(-£i, £2) A^-^l-^i, mi, . . . , m2, £2) A*-^(-^2, m2 + 1, . . . , mi - 1, 4). (3.1) 




Figure 3. The cut to be evaluated in eqn.(3.1) 



In a 1-loop amplitude with massless particles, the form of the integrals which may 
appear is well known ( see for example refs. [22,23]) and the one-loop amplitude can be 
written in the form, 

= ^ ] Cili , (3-2) 

where the coefficients Ci are rational functions of the momentum invariants and !Fn is a set 
of integral functions. In the appendices of [13] a particular choice of JF„ is given although 
many are possible [22,23]. All the cuts arise from the integral functions J^. These are 
for a massless theory, rather simple, just arising from the cuts in logarithms or the cuts 
in dilogarithms. One may see this by simple inspection of the basis of functions in [13]. 
In general one cannot reconstruct the total amplitude from the simplistic application of 
eqn. (3.1) however one may be able to reconstruct the structure of the infinities. To see 
this consider one of the integrals in the set J^n, The "one-mass" triangle which depends 
only on the momentum invariant of the massive leg, s = ^ 0, 



4? = ^(-)-'-' - (^(3 - ^ + (3.3) 



Clearly, the 1/e^ pole is closely related to the ln(— s)/e term. The latter is detectable from 
the Cutkosky rules and from it we can reconstruct the first term. This is true in general, 
by evaluating the cuts to order e° we can obtain the divergences (or rather non-cut terms 
) to order . We will use this technique to obtain the divergences present in several 
amplitudes. 

As a technical issue we shall not evaluate (3.1) but instead evaluate the the off-shell 
integral 



j r fjDp^ 1 1 

2 J (2^ ^^'^''(-^i: ms, £2)^ A'^-%-^2. + 1, . . . , mi - 1, 4)^ 



(3.4) 



cut 



whose cut in this channel is (3.1) [12]. This replacement is valid only in this channel. In 
evaluating this off-shell integral, we may substitute £f = £| = in the numerator; terms 
with l\ or £2 in the numerator do not produce a cut in this channel because the i\ or £| 
cancels a cut propagator. We emphasise that the cuts are evaluated not for a lone diagram 
at a time, but for the whole amplitude. After performing the cuts in all channels we may 
reconstruct the infinities. In practice this means evaluation the amplitude up to finite 
rational terms. 

To illustrate this consider the case of a four-point function with massless particles. 
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The set of integrals can be chosen to be fairly simple. Firstly the "scalar box integral" 




(3.5) 



together with Ij^{s,u) and I^^it^u). Then the "triangle and bubble integrals". 



/2(.) 



h{s) 




(3.6) 



Ms) 



and the corresponding functions of u and t. The function J2{s) is a linear combination 
of bubble integrals (see section 3 of ref [13]. ). In supersymmetric amplitudes the set JF4 
may be more restrictive. (Specifically J2 may be absent.) A general massless four-point 
amplitude may be written thus, 

A4^=ail^{s,t) + 02/4(3, tt) -\- a2,I^{t,u) -\-hih{s) +b2h{t) + 63/3(14) 

(3.7) 

+ Ci/2(s) +C2/2(t)+ Cshiu) +dJ2 

where the coefficients etc. are rational functions of the momentum invariants. ^ With 
this choice of JF4 it is clear which coefficients may be fixed by the cuts. For example the 
only integral function containing a ln(— s) ln{—t) term is 14(8, t) so by examining the ln(— s) 
cut we can pick out ai from the ln(— t) term in the cut. (Or alternatively the ln(— s) ln(— t) 
term if we evaluate eqn. (3.4).) Similarly by examining the other channels the are simply 
fixed. The remaining ln^(— s) contribution fixes bi and similarly the other hi. Finally the 
ln(— s) term will determine ci. So the a^, bi and Cj may be determined from the cuts with 
the only remaining ambiguity arising in the dJ2 term. For our purposes this shows that 
the only ambiguity will be in finite rational terms - and thus we may determine the infinity 
structure purely from the cuts. 

4. Example: Gravity Coupled to Scalar Matter 

In this section we will show how the results of the previous sections may be used to 
determine the UV divergences in a specific example. The example we choose is that of 

^ If we were to consider the amplitude to higher order than 0(e°) we would have to expand our basis 
to include further integral functions. 
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gravity coupled to scalar matter. The form of the UV divergences in this theory is well 
known [3] . Although, pure gravity is one-loop UV finite, in the presence of scalar matter 
infinities are generated in amplitudes which necessitate a counter-term in the action. 



203 .„„ ,,2 



320e 



{D^(t^D^^4>Y (4.1) 



Such a term is not present in the original theory and indicates that gravity coupled to scalar 
matter is non-renormalisable. We hope to arrive at the same conclusions by a consideration 
of amplitudes. The simplest amplitude effected by such a term is one with four external 
scalars ( and no external gravitons). We look at a theory with only one type of scalar. 
The intermediate states in eqn.(3.4) can be either gravitons or scalars. We consider the 
easier case of intermediate Scalars first. The tree amplitude for four scalars (all the same 
flavour) is : 

yl*-(01, 0„ 03, 04) = ^ + + i + + ^ + ^"^-^^ 

SO the product of tree amplitudes appearing in the s-channel cut is 

X (i (4(^3 -i.?^ 4(fa . eo^) + 5^ (.^ + 4(*, . 

+ 2(*r:«('<*-''*'^'') 

(4.3) 

Inserting this into eqn.(3.4) yields a variety of terms. For example the first term in the 
expansion of the above has only two propagators (l/^i and 1/^1) so may be identified 
as a tensor bubble integral and hence evaluated. In total the expansion contains boxes, 
triangles, and bubbles. These may be evaluated fairly easily using standard techniques, 
giving an expression for the logarithmic parts 

' ln(— s)- 



(47r)2-^ V 2Mt ' ' e 

+ i (3 + 2 + 2 1^) ln(-.)^ + ^'^^ ^"^-^^ ^"^-^^ + '"'^ ^"^7^ ^"^-"^ (4.4) 
4 ^ / V / 2t 2u 



+ — (-161^2 + 39 ut - 161 w^) In(-s) 1 
240^ ' ) 
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Next consider the contribution to the cut from intermediate gravitons. The tree- 
amphtudes needed are those involving two scalars and two gravitons, 



At-«(0i,2+,3+,04)=O 

.tree. . O" 3+ ^ - ^^'^^ 

where we have used a spinor hehcity convention. The spinor hehcity method for vectors is 
an exphcit reahsation of the polarisation tensors in terms of spinors 

ei^Hk;g) = ^^Ml, = i^M^, (4.6) 

^ ^ V2{q-\k+) M V '^^ V2{k+\q-) ^ ' 

where \k^) is a Weyl spinor, with plus and minus helicities, k is the on-shell momentum 
of the vector and q is an arbitrary reference momentum satisfying = 0, /c ■ g 7^ 0. The 
spinor helicity method for gravitons [25,26] is related to that for vectors [27] by 

where are the graviton helicity polarisations and £^ are the vector helicity polarisations 
defined by Xu, Zhang and Chang. We use the notation for spinor inner products |/c^) = 
(12) and {k'l\k2) = [12]. The use of spinor helicity techniques has proved extremely useful 
in QCD calculation and we will take advantage of the benefits here also. All states are 
taken to be outgoing and may have plus or minus helicity. 

Equation (4.5) implies that the only contribution comes from cuts with gravitons of 
differing helicity across the cut. The contribution to the cut from these will be 

(01 , <^2 , , (-t ) A'^^^ {<f>s,<t>4,i2Jt) 

(l£i)"(2£i)" s{k^-h) (4^2)' (3 £2)" s{h-h,) 

16((14)(^1^2)(^2 2)(21))2 {k^-l2) ((3£2)(^2^l)(^l4)(43))2 (£i./c3) 

{tr-{hhhh)Y 

16 s2 (2£i • kif{2ii ■ k2){2h ■ k^){2h ■ k^) 

(4.8) 

Where [ah] {be) [c d] ■ ■ ■ {m a) ~ tr_(^^^- • -fn) and tr±(p) = tr(-^i^^)p) (We must also 
include the contribution with the other choice of helicities on the internal legs: 

A{cj>^,ct>2,it,^2)Mh,<l>^,4,^l) (4.9) 
This is equivalent to setting £2 — > —£i and £1 — > —£2 ) 
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If we combine these two contributions and carry out the integrations (which a rather 
more comphcated than those for the scalar case) we find 



K 



(47r) 



2-e 



/{u^ -ut + r) ^ ^ ^„ 1 , , ^ ^ ^ 
^ H-sf - ^ + + ^ ) ^^(-^) 

^ In(-t) In(-s) ^ In(-u) In(-s) 



4s 4s 
If we add this to the result for intermediate scalars we obtain 
K^rr + 2 uH + 2ut^ + Z uH^ + t^) In(-s) 1 



(4.10) 



(47r)2-«\ 2ut e 

+ ^{3u^ + 2ut + 3t^) rr ln{-sf 



(^4 + s4) In(-s) In(-t) (s^ + «4) ln(-«) ^"^'^^^ 



+ ^ — ^ — - + 

+ 



2ts 2 us 

1 {-l%ZuH-43,uH'^ -lQ2,ut^)\a{-s) 



240 ut 

From this we can write down an expression with the correct cuts in all channels 

8(s^ + t^)h{s, t) + 8(s^ + u'^)h{s, u) + 8('u^ + t'^)h{t, u) 

- 8s(3s2 + + u^)h{s) - 8t{3t'^ + s'^ + u'^)h{t) - Su{3u^ + + s'^)h{u) 
2(163^2 + 163^2 + 43tu) ^ , , 2(163^2 + 1635^ + 43^5) ^ , , 

+ ^^^'^ + 15 ^''^^^ 

2(163s2 + 163*2 + 43ts) ^ , , 
+ h{u) 

(4.12) 

This expression has the correct cuts in all channels and contains the correct infinite pieces. 
It is not the correct full result. As discussed in the previous section we expect the full answer 
to differ from the above by finite, non-logarithmic rational polynomials in the momentum 
invariants. However the above expression does contain the correct 1/e divergences. 
Contained in the above we expect a IR infinity of the form 



irr^' ((-.)^- + (-t)^- + (-^)^-) . . .^ 

A (01,02,03,04) 



2(47r)2-^ e2 

( 



rr«4 ((-s)i-^ + {-tf-' + {-uf-') (s'^ + e ^ s2+^2 + u'' 



(4.13) 



2(47r)2-^ e2 \ 2u 2t 2s 

By examination we can identify this term plus an additional infinity, 

-||('^+«^+"^) (4.14) 
12 



Since the IR divergences are accounted for, this must be a UV infinity. By inspection, we 
can see this corresponds to the counter term in the Lagrangian in eqn. (4.1) - with the 
correct coefficient. 

Before finishing this section we can consider the case where gravity is couple to a set 
of N scalar 0^. Using the techniques above we can identify the counter Lagrangian to be 

Y, ( {^^^^) (s^<t>iS''-l>i)(aAia"<t>i) (4.15) 

It is possible to examine the counter terms for a variety of types of matter coupled to 
gravity. Further examples will be given elsewhere [28] . 



5. The Cut Calculation of Al-loop(i+ 2+, 3+, 4+) with a massive scalar. 

Although a naive interpretation the Cutkosky rules suggests that one may only use 
them to evaluate the amplitudes up to finite rational polynomials in the momentum in- 
variants it may be possible to use the rules to evaluate these also. The idea is fairly simple. 
In dimensional regularisation, if we can evaluate the cuts to order then we can recon- 
struct the rational polynomial terms to order e"^"^. The difficulty lies in evaluating the 
cuts consistently to all orders. That this is possible, has been demonstrated by Bern and 
Morgan in ref. [14] where it was shown how to, with care, evaluate the cuts to higher order 
in e. Previously we have used the on-shell tree amplitudes with intermediate legs having 
momenta in D = 4. However the loop momentum integral has momenta in Z) = 4 — 2e. 
This involves an error, which although not contributing to finite order in ln(— s) [13], gives 
an error in terms eln(— s) which feeds down to finite polynomial terms. To correctly carry 
out the Z) = 4 — 2e loop momentum integrals it is convenient to split the Z)-dimensional 
momentum into a four dimensional and a — 2e momentum, whence 

pI^pI-I^^ (5.1) 

and the integration splits up as 

d^p d^p d-^'ii (5.2) 

This prescription is well described in ref. [29] . If calculating the amplitude for a massless 
scalar circulating in the loop the prescription is now clear: one uses on-shell four dimen- 
sional amplitudes but where the scalar now has a mass //^ and one integrates according to 
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eqn. (5.2). Hopefully an example will make this clear. The first example we consider will 
be the contribution to four graviton scattering when the external gravitons all have the 
same helicity from a massless scalar circulating in the loop, 

4Vair(1^2+,3+,4+) (5.3) 

^ This can be found in ref. [10] where it was calculated using string-based diagrammatic 
rules. The answer is a finite rational polynomial and hence we might expect it not to be 
calculable using cuts. This is consistent with the fact that the tree amplitudes necessary 
to evaluate the cuts A^^'^'^{1~^ , 2+, cp, (j)) vanish for a massless scalar. However if we consider 
the tree amplitude for a scalar with A mass ji^ we find, 

A,{1\ 2+, 3+, 4^) = _^^ (^') [2 3]' / 1 1 \ ^5 4^ 

When we calculate the cut in the, for example, Si2 channel the two factors yield two box 
integrals with ordering of legs 1234 and 1243. The 1234 ordering will contribute (after 
inserting the factors of jj?) and 1/(^1 " A*')) 

[12]2[3 4]2 /• d^p 



16(12)2(34)2 y (27r)4-2e 



X 



(p2 - jj?){{p - /Ci)2 - IJ?){{p -ki- /C2)2 - /x2)((p - /Cl - /C2 - ks)'^ - jJ?) ' 

(5.5) 

The integrals over — 2e loop momenta are evaluated using eqn. (22) of ref. [29]. 

and (after Feynman parameterising the J df^p integral ) 

- diJ,\i?y^-'- 2-2 = ""f^ ~ f {-saias - ta2a^f-' (5.7) 

We can then reconstruct the cut in the s-channel 

iK^ [12]2[3 4]2 7re(3-e) 6(1-6) (2-e) /,,,=i2-2e , .^=12-2; 
16(47r)2-^ (12)2(34)2 sin(7re) r(l - e)r(l + e 

where the definition of the D-dimensional box integral is 



ID ^ _ ^/2) f da, S{Ea^ - 1)-, ] 



(5.9) 



t Due to supersymmetric Ward identities [31] this amplitude due to massless scalars circulating is equal 
to that for gravitons circulating. 
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By symmetry we can deduce that the following object has the correct all- order in e cuts 
in all three channels and hence must be the correct amplitude, 

iK^ [12]- [3 1]- C 7^(3 - c(l - r)(2 - c) |^^j^^i2-2e ,TD=12-2e ,rD=12-2e 

16 (12)^(34)M47r)2-^r(l-6) sin(7r6) T{1 + e) [''^'^ ^'^^^^ 

(5.10) 

Since this is now the correct answer to all order in e A we can evaluate it merely at order 
e°. Thus we need /'D=i2-2e order 

tD=12-2€ _T^( r, , \ ,„ -IM' „ + \2-e 



= r(-2 + ^) J SC^tti - l)(-saia3 - ta2a4y 
= Ljdai 5{Y,ai - l)(-saia3 - ta^a^)"" + 0(e°) (5.11) 



5040e ^ ' 

and thus the amplitude is 

[1 2f [3 Af 6 4s2 + At^ + Au'^ + st + tu + su 



(12)' (34)' 16(47r)2 2520 

/ St \2(s2 + t2 + ^2) 



(5.12) 

It") 

(47r)2 V(12) (2 3) (34) (41)7 3840 
This reproduces the result of ref. [10] for a real scalar. (For a complex scalar multiply by 
a factor of 2). and is consistent with [30]. thus we have thus demonstrated, as was done 
for QCD in [14] how the cuts can reproduce the finite rational polynomials in amplitudes. 

After demonstrating the validity of the method let us calculate the amplitude with 
a massive scalar of mass M in the loop. This has been calculated in the M oo limit 
previously [16]. This will involve replacing ^"^ by + in the initial cut expression. 
Following the analysis to eqn. (5.5) /x^ is replaced by {iJ? + M^)^ in the numerator and jj? 
by jj? + in the numerators. The various terms can then be evaluated using 

/j-oo 2n 



{saias + to2a4 - M^)"-^-^ (5.13) 



(— l)"'7r(n — e — 1 1 , ,^9n„ ^ 9 

sin(7re) 

1 rD=2n-2e+A 



sin(7re) r(2 + e - n) 

Putting the pieces together we may obtain 

iK^ [12]' [3 4]' ire / -e)(2- e)(3 - eU^=i'-'^ 

16(47r)2-(12)'(34)'sin(7r6)r(l-6)r(l + 6)lv '^^^'^^ 

+ 4M'6(1 - 6)(2 - 6)/f2=f -'^ + 6M^6(1 - e)/f2=|-'^ (5.14) 
+ AM^el^^f-^^ - M^I^^m'^' + {1243} + {1324} 
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which will be valid to all orders in e. To evaluate to e°, we will only need the infinite parts 
of the D = 12,D = 10 D = 8 and D = 6 integrals and the finite part of the Z) = 4 integral. 
(The D = Q integral will drop out to order e° since the box is both UV and IR finite in 
D = 6.) We can evaluate these objects to 0(e°), using 

<}^0^I^=l^-^^A = J da, 5{j:a, - l){sa,as + ta^a, - M^f'^ 

_ + st + 2t^ _ M^{s + t) 
~ 2520 60 ^ ~6" 



r{l + e) ^^^^^~^' " / ^^^""^ ~ l)(soia3 + to2a4 - M^)"^ = 1/6 
which will yield for the amplitude (for a real scalar), 



(5.15) 



/ st V / (s'^ + 1'^ + u'^) 1, . 



16(47r)2 V(12) (2 3) (3 4) (41)7 V 240 2 

-'^^ 1234 ^ + -'1243 +-'1324 ^ 



(5.16) 



This amplitude contains logarithms and dilogarithms from Ii2si- These terms could be 
correctly obtained from the cuts naively but we also have the additional polynomial terms 
to give the full correct result. We might also expect that the amplitude should contain 
terms of the form hi{m'^), which would not have been recovered in the above analysis. 
However, we can deduce that no such terms appear by applying the arguments outlined 
in ref. [14]: Such terms only appear with specific divergent contributions; since we know 
that the amplitude is UV and IR finite we know that no ln{m^) terms will be found. The 
large M expansion, on /^^4 is 



J da, S(^ai - l){saia3 + ^204 - M^) 



2\-2-e 



/ 



da, 5(Ea. - 1) (M^) ( 1 - (2 + e) ^ + (2 + e) (3 + e) ^ 



- (2 + 6)(3 + e)(4 + 6)^ + (2 + 6)(3 + 6)(4 + 6)(5 + 6)^ (5.17) 

- (2 + 6)(3 + 6)(4 + 6)(5 + 6)(6 + 6)^2!^ • • • ) 
1 s + t _^ 28"^ + st + 2t'^ 3s^ + st^ + s'^t + 2,t^ 



6M4 60M6 840M8 7560Mio 
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where S = saia^ + ta2a4. From this we obtain the large mass expansion for the amphtude 

~ 16(47r)2 V (1 2) (2 3) (3 4) (4 1) ) \bMM^ 15120M4 15840M6 

(5.18) 

This amphtude tends to zero as M — oo as expected with the 1/M^ term having the same 
coefficient as calculated previously [16]. 



6. Conclusions 

In perturbative gravity and gauge theories the calculation of amplitudes is a painful 
process where huge intermediate expression eventually collapse to relatively small final 
answers. This computational complexity will remain a challenge despite considerable im- 
provement in both techniques and algebraic computing facilities. Any technique which 
attempts to avoid this computational explosion must be explored carefully. Unitarity in 
the form of the Cutkosky rules is one such technique. By sewing together tree amplitudes 
upon which much simplifications have already been performed one can build upon previous 
calculations and avoid part of the algebraic complexity. Unfortunately, the full amplitude 
is not calculated by the Cutkosky rules but the rules are a powerful constraint upon the 
form of amplitudes. In some cases unitarity is enough to fix the the amplitudes completely, 
however naively applying the Cutkosky rules in the general case leaves ambiguities. These 
ambiguities however can be expressed in such a form that they only effect the finite terms 
polynomial in momentum invariants. Thus, in principle, one can calculate the infinities 
using the Cutkosky rules. More sophisticated use of the Cutkosky rules [14] can also be 
used to determine the entire structure of a loop amplitude. In quantum gravity it is most 
often the infinities in amplitudes which are of more interest rather than the finite terms 
(The case in QCD is often the other way aroimd.). Here we have shown in several cases 
how the Cutkosky rules can be practically used to evaluate infinities. In any massless 
theory there are both IR and UV infinities. Although the cuts do not distinguish between 
these sources, the general structure of the IR infinities can be determined [16]. With this 
information we can isolate the UV divergences. We have explicitly shown how the UV 
infinities in scalar coupled matter appear in an amplitude using this technique and are in 
agreement with the known results which were obtained by a very different route. Other 
infinities previously unknown may be obtained in this manner [28]. 

We would like to thank Zvi Bern and Andrew Morgan for useful discussions regarding 
the use of Cutkosky rules for massive particles. Christian Schubert for discussions relating 
to the world-line formalism for gravity and Lance Dixon for discusions both regarding the 
possibility of disentangling the UV and IR infinities within the amplitudes and for access 
to his calculations regarding massive scalar amplitudes. This research was supported by 
PPARC. 
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